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The Fock-Darwin States of Dirac Electrons in Graphene-based Artificial Atoms 
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We have investigated the Fock-Darwin states of the massless chiral fermions confined in a graphitic 
paraboUc quantum dot. In the light of the Klein tunneling, we have analyzed the condition for 
confinement of the Dirac fermions in a cylindrically-symmetric potential. New features of the energy 
levels of the Dirac electrons as compared to the conventional electronic systems are dicussed. We 
have also evaluated the dipole-allowed transitions in the energy levels of the dots. We propose 
that in the high magnetic field limit, the band parameters can be accurately determined from the 
dipole-allowed transitions. 
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Quantum dots (QDs), or the 'artificial atoms' [l[ are 
one of the most intensely studied systems in condensed 
matter physics where the fundamental effects related to 
various quantum phenomena in confined geometries can 
be studied but with the unique advantage that the na- 
ture of the confinement and the electron density can be 
tuned externally. However, much of the interest on this 
system derives from its enormous potentials for applica- 
tions, ranging from novel lasers to quantum information 
processing. While the majority of the QD systems inves- 
tigated are based on the semiconductor heterostructures, 
in recent years, quantum dots created in the carbon nan- 
otubes have been reported in the literature where the 
'atomic' properties were clearly elucidated and its im- 
portance in technological applications was also demon- 
strated . Conductance properties of ultrathin graphitic 
QDs [1] have also been reported recently. It is now 
well recognized that the low-energy dynamics of the two- 
dimensional electrons in graphene is governed by the 
Dirac- Weyl equation, and the charge carriers behave as 
massless chiral fermions Q . In this situation, confine- 
ment of electrons becomes quite a challenging task, due to 
the so-called Klein's paradox [3| . This problem has been 
dealt with in the case of a one-dimensional (ID) wire in 
zero and in finite Q magnetic fields. In this letter, we 
report on the electronic properties of the parabolic QDs 
in graphene, in particular, we present the energy levels 
as a function of the magnetic field (Fock-Darwin states 
[H) and the associated dipole-allowed optical transitions 
in this system. We propose that the optical spectroscopy 
of the graphene QD in the high-field limit could provide 
an accurate means of determining the band parameters 
of graphene. 

The Hamiltonian of a single electron in graphene with 
a cylindrically symmetric confinement potential is 

H = Ho + Hi = J(a7r) + V^(r), (1) 

where a are the Pauli matrices, tt = p + |v4, A = 
■j(— y, x) is the vector potential corresponding to the 
magnetic field B in the z-dircction orthogonal to the 



graphene plane, and 7 = -y/S 070/2 is the band param- 
eter. Here a = 0.246 nm is the lattice constant and 
7o (meV) is the transfer integral between the nearest- 
neighbor carbon atoms [lo| . 

At first we analyze the properties of the graphene sys- 
tem in the absence of a magnetic field to find the condi- 
tion for confinement of an electron in the potential V{r). 
Due to the Klein tunneling the electrons in graphene can 
not be localized by a confinement potential, since for any 
potential there will be the electron states with negative 
energy (the hole states) which would provide the escape 
channel for the electron inside the potential well. We can 
then discuss only the quasilocalized states or trapping of 
the electron by the confinement potential. This problem 
has been treated for the quasi-lD graphene system [s, 11|, 
where it was shown that the transverse momentum in ID 
introduces the classically forbidden regions, which helps 
in trapping the electron. The width of the quasilocalized 
level is determined by the tunneling through the classi- 
cally forbidden regions. For the zero transverse momen- 
tum, the tunneling barriers disappear and there are no 
trapped states. In our case, we have a cylindrically sym- 
metric confinement potential with the effective transverse 
momentum m/r, where m is the electron angular momen- 
tum. Therefore, for m 7^ we expect the trapping of an 
electron by a cylindrically symmetric QD. In terms of the 
two-component wave function (xi(''')e'^™~"'^-'^, X2(?')e""'^) 
the Schrodingcr equation corresponding to the Hamilto- 
nian ([1]), is 



dx2 ™ 

V{r)xi - tj-: «7— X2 = Exi 

ar r 

^n \ . dxi TO - 1 

ar r 



(2) 
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There are no analytical solution to these equations, and 
therefore, we first present below a semiclassical analysis. 

Semiclassical analysis: At a large to we seek a solution 
of Eqs. ([SJ-jS]) in the form e*^'', which gives 
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The classical turning points can be found from the condi- 
tion g = 0, and the classical region is y(r)| > ^\m\/r. 
If ro is the solution of the equation E — V{r) = then we 
can find the classically forbidcn region as (ro — Ar) < r < 
(ro + Ar) , where Ar ~ m/Fro, F ~ "f^^dV{ro)/dr, and 
we assumed that F ^ m/r'^. If the electron is trapped 
in the dot, i.e., at r < 7'o — Ar, then the escape rate or 
the width of the quasilocalized levels is determined by 
the tunneling through the classically forbidden region. 



T = exp - 



ro + Ar 



ro — Ar 



\q{r)\dr 
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Therefore, in order to trap the electron we need a large 
m and a small F , i.e. a smooth confinement potential. 
For a potential V = {u/n)r^\ Eq. ([5]) takes the form 



exp 
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exp 



n(i;/e„)"+i/'' 



(6) 



where e„ — (7"w/ri,)^/("+^^. Equation ([6]) gives the 
upper limit on the energy of the quasilocalized levels 
at a given m, i.e. Ejtn < m^"/^""'""'^^. Based on the 
semiclassical expression we can also find the interlevel 
separation of the quasilocalized levels at large energies, 
AEn = Q;e„(i?/e„)~^/", where a ~ 1. We then estimate 
the number of quasilocalized levels, Nn,m for a given an- 
gular momentum, m, and a given potential profile from 



dE 
A^ 



1 



(7) 



This estimation is valid for a large m. For a small m we 
need to solve the system of equations Q-© numerically 
to find the properties of the quasilocalized states. 

Quasilocalized states - numerical solutions: To extract 
the information about the width of the quasibound levels 
we need to impose special boundary conditions far from 
the QD. This condition means that far from the origin, 
r 3> foj the solution should be an outgoing wave, i.e., the 
propogation away from the QD. From Eqs. (HI)-®, it is 
clear that for r ^ tq the outgoing solution has the form 



Xiir) = -X2{r) = C exp 
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V{r')dr' 



(8) 



where C is a constant. Equation ([5]) is the boundary 
condition for the system (H))-© at large distances. Since 
at r = the solution should be non-divergent, another 
boundary condition is = 0) = X2{r = 0) = 0. A 

solution with these boundary conditions exists only for 
a complex energy, E. The imaginary part of the energy 
determines the width of the quasilocalized level. 

We have solved Eqs. (H))-® numerically for a poten- 
tial of the form V{r) ~ with different values of 
the exponent, rt, and for different values of the angular 
momentum m. In the dimensionless units, i.e., for the 
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FIG. 1: The real and imaginary parts of the energy spectra 
of an electron in a QD with a confinement potential V{r) = 
(u/?i)r", shown for various values of the exponent n and the 
angular momentum m: (a) n = 2, m = 2 (circles), and n = 2, 
m = 10 (dots); (b) n = 4, m — 2 (circles), and n — 4, m = 10 
(dots). The results for n = 2 and m — 1 are shown as inset. 
The energy is in units of e„. 



units of length and the energy, a„ = (717/1/)^/'""''^' and 
= (7"M/ri,)-'^/("+^^ respectively, the system does not 
contain any information about the interaction strength, 
u. The results in Fig. 1 show a bunch of closely spaced 
levels with a large imaginary part of the energies, Im[E). 
For these states Im{E) is comparable to the interlevel 
spacing. These states are the delocalized contunuum 
states. At the background of the continuum spectra we 
also see the levels with a small imaginary energy, i.e. a 
small width of the level. These states are the quasilocal- 
ized states of the QD. The manifestation of such states 
can be seen already at m = 2 (circles) both for n = 2 
[Fig. 1(a)] and n = 4 [Fig. 1(b)]. The strength of the 
localization can be characterized in terms of the ratio 
of the Im{E) to the interlevel spacing. For m = 2 this 
ratio is 50. With an increase of the magnetic moment 
the quasilocalized states become well developed and at 
m = 10 we clearly see the states with very low Im{E). 
The ratio of Im{E) to the interlevel spacing for these 
states is about 800. With an increase of the energy the 
states are less localized, i.e. Im{E) increases. This is 
consistent with Eq. Note that for all values of the 
exponent n there are no localized states at m ~ 1 [inset in 
Fig. 1(a)]. All the states at to = have very large Im{E). 
There are also no quasilocalized states at m = 0. The 
reason for delocalization of the electron at to = and 1 is 
that the effective transverse momentum for either the xi 
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component (at m = 0) or the xi component (at m = 1) 
is zero. In the following, we analyze the magnetic field 
effects on the electronic states of the QDs. 

Magnetic Field - semiclassical analysis: In a magnetic 
field the system of equations (HI)-® has an additional 
non-diagonal term proportional to the magnetic field. In 
the dimcnsionlcss units, i.e., for the units of length a„ 
and the energy e„. the system is characterized by only 
one parameter, b = {eB /2c){n'-f /u)''^^^"''^^\ In the semi- 
classical approximation, the effective transverse momen- 
tum is (m/r + br) and the Eq. Q becomes 



T = exp 



7r(m + 5£;2/")2 



(9) 



where E = i?/e„. The effect of the magnetic field is 
different for the states with a positive or a negative m 
(the sign of m depends on the direction of a magnetic 
field). For a positive m the application of a magnetic 
field increases the effective transverse momentum and 
suppresses the tunneling from the QD. For a negative 
TO, the magnetic field decreases the tranverse momen- 
tum. Therefore the state becomes less localized. If we 
increase the magnetic field even further then at some 
point, b = m/E^^", the level becomes delocalized, and at 
a even larger B the level again becomes localized. Now 
the trapping will be due to the magnetic field. Therefore 
for a negative m, the magnetic field induces a localization- 
delocalization-localization transition. 

The number of the quasilocalized states in a weak mag- 
netic field is estimated to be Nn.m '--^ [to -I- 6|TOp/^"+^^] . 
This number with a positive m increases with an increas- 
ing magnetic field, while that for a negative m decreases 
with the magnetic field upto a certain value of B and 
then increases. The total number of states with posi- 
tive and negative angular momenta, Nn^m + Nn^-m ~ 
771^ + 5^|to|''^/'^""''^) always increases with an increasing B. 
From this bahavior we expect the following effect: We 
assume that the QD is occupied by electrons upto a cer- 
tain energy, i.e., the states with both positive and nega- 
tive angular momenta are occupied and the net angular 
momentum of the dot is zero. We then apply a mag- 
netic field and the states with positive m, becomes more 
localized while the electrons from the states with nega- 
tive TO will escape from the QD. Finally, the electrons in 
the QD will have a net positive angular momentum and 
correspondingly a net magnetic moment. 

Magnetic Field - numerical results: To study the de- 
pendence of the quasilocalized spectra on the magnetic 
field we introduce the wavefunctions of the Hamiltonian 
Ho, i.e. without a confinement potential |I2l |. as the ba- 
sis functions. To eliminate any escape of the electron 
from the QD we consider only the basis functions with 
the positive energy, 



,?n — 1 (^) 
i4>n,m{x) 



where N ~ n 

index, Cn=q = 
and 

1 



f ^(|to| + m) is the Landau Level (LL) 
f and Cjv#o = l/\/2, sgn{N = 0) = 0, 



V2 V ("+ \m\y- 



e'^/\l™l/24™l(a;)e™'' (11) 



is the Landau wavcfunction. Here ill" '(a;) is the associ- 
ated Lagucrre polynomial, x = r'^/a' is a dimcnsionlcss 
distance. Here a' is the characteristic length of the sys- 
tem. Without the confinement a' should be equal to the 
mangctic length / = he/ eB. In the presence of the 
confinement, the Hamiltonian suggests a natural unit of 
length {'-^ /u)"^/^ Q and a natural unit of energy {'^'^uY/^. 
This length characterizes the size of a parabolic dot in 
graphene. Therefore, 



2 



1 

/2 



2/3 




(10) 



B (Tesia) 



FIG. 2: Fock-Darwin spectrum of the Dirac quantum dots, 
plotted for the confinement potential strength m = 0.1 (in 
dimensionless units). The numbers in the parentheses corre- 
spond to the two quantum numbers n and m. Results for 
u = are given as an inset. 

In our numerical calculations, we choose the band pa- 
rameter to be 7 646 meV-nm for 70 = 3.03 eV jl3| . 
The low-lying energy states of the graphene QD are 
shown in Fig. 2. In the absence of a confinement poten- 
tial, the Dirac spectrum scales as ^/2^^fN /I [shown as 
inset]. In the Fock-Darwin spectrum for a conventional 
electron dot, the energy levels are degenerate and equally 
spaced at i? = [H . The two-dimensional parabolic con- 
finement considered here shows two outstanding features 
in contrast to the Fock-Darwin spectra at B = 0. The 
first is the lifting of the degeneracy and the other is the 
unequal seperation among the energy levels. Figure 2 
shows the field-dependent energy spectrum for u — 0.1. 
The energy difference between the lowest two levels at 
_B = is about (7^^)^/'^. At a low magnetic field, the 
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magnetic length I is larger than or comparable to the size 
of the confinement (7/1*)^/'^ and there is a hybridization 
of the LLs with the levels arising from the spatial con- 
finement. In the high magnetic field limit I <C {■y/uy^^, 
the Landau-type levels prevail, as expected. The Fock- 
Darwin spectra for conventional quantum dots have been 
determined earlier by the transport spectroscopy [14|. 



Similar studies for the graphene QDs would be very im- 
portant to explore the energy levels and the nature of 
confinement for Dirac fermions in a graphene QD. 
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FIG. 3: The dipole-allowed transitions in the Fock-Darwin 
spectrum of graphene QDs for it = 0.1. Inset: the case oiu = 
0. The thickness of the lines is proportional to the calculated 
intensity. From the bottom to the top, the relative intensities 
are about 1.0, 0.1, 0.02, respectively. 

Figure 3 shows the dipole-allowed optical absorption 
spectra [l|, [I^ for u = and u = 0.1. Without any 
confinement there is only the (0, 0) (0, 1) transition 
(shown as inset in Fig. 3). The additional transitions, 
(0, 0) (1, 1) and (0, 0) (2, 1) are due to the presence 
of the parabolic confinement. For the Dirac fermions in 
a QD, the lowest dipole-allowed transition is 
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For 5 = the corresponding energy is sa [j'^uf-l'^ . As the 
magnetic field increases, AiJ approaches the cyclotron 
energy ^/2(e7B/ftc), where 7 can therefore be 

uniquely determined experimentally. Conventionally, in 
the nearest-neighbor tight-binding model, 70 is obtained 
by fitting the ab initio calculation and the experimen- 
tal data [l3|- In a GaAs quantum dot, the magnetic- 
field-dependent far-infrared absorption experiments have 
established the e nerg y relation /S.E± = Ml ± ^hiUc to 
a great accuracy [l5|. Similarly, for the masslcss chiral 
fermions in a graphene QD, we expect that the band pa- 



rameter 7 can also be determined quite accurately by the 
optical absorption experiments in the high-field limit. 

Although the high magnetic field results are the major 
focus of this paper as in this case the localization of the 
electron in a QD at all values of the m is provided by 
the magnetic field, at B=0 the chiral nature of the states 
prevent the electrons from being confined in a QD. This 
clearly indicates that the nature of the energy states and 
the optical spectra at a very small B arc still important 
open questions. 
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